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HIGHLIGHTS 


•  This  is  the  first  attempt  on  studying  pyrolysis  kinetics  of  Polysiphonia  elongata. 

•  Kinetic  parameters  were  determined  from  TG  analysis. 

•  Pyrolysis  process  was  modeled  using  obtained  kinetic  parameters. 

•  Model  was  in  a  good  agreement  with  experimental  data. 

•  Results  can  be  helpful  for  design  of  pyrolytic  conversion  systems. 
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The  pyrolysis  characteristics  and  kinetics  of  Polysiphonia  elongata  were  investigated  using  a  thermogravi- 
metric  analyzer.  The  main  decomposition  of  samples  occurred  between  225  °C  and  485  °C  at  heating 
rates  of  5-40  °C/min;  owing  to  release  of  78-82%  of  total  volatiles.  The  heating  rate  effected  pyrolysis 
characteristics  such  as  maximum  devolatilization  rate  and  decomposition  temperature.  However,  total 
volatile  matter  yield  was  not  significantly  affected  by  heating  rate.  The  activation  energy  of  pyrolysis 
reaction  was  calculated  by  model  free  Friedman  and  Kissenger-Akahira-Sunose  methods  and  mean  val¬ 
ues  were  1 16.23  kj/mol  and  126.48  kj/mol,  respectively.  A  variance  in  the  activation  energy  with  the  pro¬ 
ceeding  conversions  was  observed  for  the  models  applied,  which  shows  that  the  pyrolysis  process  was 
composed  of  multi-step  kinetics.  The  Coats-Redfern  method  was  used  to  determine  pre-exponential 
factor  and  reaction  order.  The  obtained  parameters  were  used  in  simulation  of  pyrolysis  process  and 
results  were  in  a  good  agreement  with  experimental  data. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Fossil  fuels  are  the  main  energy  sources.  However,  they  are  lim¬ 
ited  and  the  depletion  of  these  energy  sources  is  expected  by  the 
middle  of  this  century  (Kraan,  2013).  Furthermore,  dependence 
on  fossil  fuels  has  largely  increased  the  emissions  of  greenhouse 
gases,  particulate  matters,  and  other  pollutants,  which  results  in 
many  environmental  problems  such  as  global  warming  and  atmo¬ 
spheric  pollutions.  The  future  supply  security  of  energy  sources  is 
another  crucial  issue.  Thus,  there  has  been  an  increasing  interest  in 
alternative  and  sustainable  energy  sources.  Biomass  is  one  of  the 
most  outstanding  options  as  an  alternative  fuel  because  it  is  clean, 
renewable  and  fast  growing  (Kim  et  al.,  2013). 
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Macro  algae  have  been  proposed  to  have  great  potential  as  bio¬ 
mass  and  they  are  primarily  composed  of  polysaccharides  that  can 
be  converted  into  fuels  such  as  bio-alcohols  or  bio-oils  by  fermen¬ 
tation  or  pyrolysis,  respectively  (Kim  et  al.,  2012,  2013).  This  type 
of  biomass  is  attracting  the  attention  of  researchers  for  its  potential 
as  an  environmental-friendly  and  economically  sustainable 
resource.  The  main  advantages  include  easy  of  large-scale  produc¬ 
tion,  rapid  growth,  marine  environment  protective  effects.  Addi¬ 
tionally,  their  production  does  not  require  land  or  fresh  water  (Li 
et  al.,  2011;  Kraan,  2013;  Kim  et  al.,  2014a, b). 

The  pyrolysis  is  the  thermal  degradation  of  materials  in  the 
absence  of  oxygen  (Kim  et  al.,  2014a, b)  and  it  is  an  attractive 
way  to  use  the  energy  contained  in  the  biomass.  Solid  biomass 
and  wastes,  which  are  very  difficult  and  costly  to  manage,  can  be 
readily  converted  into  liquid,  gas  and  charcoal  products  by  the 
pyrolysis  process. 
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Almost  all  types  of  biomass  can  be  used  as  feed  stocks  for  ther¬ 
mochemical  conversion.  However,  thermochemical  behavior  of 
various  aquatic  biomasses  is  largely  different  not  only  from  ligno- 
cellulosic  biomass  but  also  from  each  other,  and  the  bio-oil 
produced  also  varies  largely  as  a  result  of  the  differences  in 
components  and  operation  conditions  (Kim  et  al.,  2013).  Thus,  a 
comprehensive  knowledge  of  the  pyrolysis  mechanism  and  kinetic 
analysis  of  various  marine  biomasses  is  required  (Wu  et  al.,  2014). 

In  addition  to  this,  both  the  development  of  the  pyrolytic  pro¬ 
cess  and  reactor  design  require  complete  clarification  of  the  pyro¬ 
lytic  mechanism.  Therefore,  pyrolysis  kinetics  and  thermal 
decomposition  mechanisms  for  macroalgae  should  be  extensively 
studied  (Li  et  al.,  2011,  2012). 

In  their  recent  study,  Wang  et  al.  (2013)  studied  the  pyrolysis  of 
macroalgaes,  Enteromorpha  clathrata  and  Sargassum  natans  and 
compared  the  composition  of  biooil  products.  Yanik  et  al.  (2013) 
published  a  study  on  algal  biomass  pyrolysis  and  reported 
pyrolysis  products  of  Laminaria  digitata  and  Fucus  serratus. 
Ferrera-Lorenzo  et  al.  (2014a, b)  investigated  the  pyrolysis  charac¬ 
teristics  of  the  algae  meal  which  was  generated  by  the  industrial 
production  of  Agar-Agar.  Kim  et  al.  (2012,  2013)  investigated  the 
alga  Sagarssum  sp.  biomass  and  Saccharina  japonica  pyrolysis  and 
reported  kinetic  parameters  using  reaction  models.  Wu  et  al. 
(2014)  studied  the  pyrolysis  of  three  different  aquatic  biomass  spe¬ 
cies  (microalgae,  macroalgae  and  duckweed)  and  calculated  kinetic 
parameters  by  using  the  thermogravimetric  analyzer. 

Turkey  is  surrounded  with  seas;  therefore  there  is  a  huge  aqua¬ 
tic  biomass  potential.  There  is  no  evidence  in  literature  of  the  iden¬ 
tification  of  kinetic  parameters  in  the  inert  atmosphere  for 
Polysiphonia  elongata.  In  this  work,  the  pyrolysis  behavior  and 
kinetic  parameters  including  apparent  activation  energy,  pre-expo¬ 
nential  factor,  and  reaction  order  of  P.  elongata  (PE),  were  investi¬ 
gated  using  TGA.  Kinetic  constants  for  the  P.  elongata  pyrolysis  in 
the  inert  atmosphere  obtained  with  different  model-free  methods 
were  not  available  in  the  literature  and  the  results  of  this  study  can 
provide  useful  knowledge  to  pyrolysis  researchers  and  engineers. 
Furthermore,  these  data  are  expected  to  help  design  and  scale 
up  the  thermochemical  conversion  processes  for  macroalgae 
P.  elongata. 


2.  Methods 

2.7.  P.  elongata  (PE) 

Alga  P.  elongata  (PE)  was  collected  from  Samsun  Harbor,  Turkey. 
The  samples  were  cleaned  in  fresh  water  followed  by  distilled 
water,  dried  in  oven  at  60  °C  overnight.  The  dried  samples  were 
pulverized  and  sieved  for  63-125  pm  of  particle  size  and  used 
for  experiments.  PE  was  characterized  in  terms  of  proximate  and 
ultimate  analysis,  according  to  the  ASTM  standards  (E871, 
D1 102-84)  by  using  a  programmable  ash  oven.  Leco  type  analyzer 
CHNS-932  was  used  for  the  elemental  analysis  and  oxygen  bomb 
calorimeter  was  used  for  the  determination  of  high  heating  value 
of  PE. 


2.2.  Thermogravimetric  analysis 

The  TGA  of  PE  was  conducted  using  a  Simultaneous  Differential 
Thermogravimetric  Analyzer  which  contains  a  heat-flux  type  DTA 
and  a  TGA  (Shimadzu,  DTG-60,  Japan;  with  a  precision  of  temper¬ 
ature  measurement  ±0.1  K,  DTA  sensitivity  ±0.1  pV  and  microbal¬ 
ance  sensitivity  ±0.1  pg).  The  instrument  provided  continuous 
recording  of  TGA  and  derivative  thermogravimetric  (DTG)  analysis 
curves  in  terms  of  weight  loss  percent  per  second.  In  each  experi¬ 
ment,  10  mg  of  sample  was  loaded  into  the  platinum  crucible  of 


the  thermal  analyzer.  To  maintain  the  inert  condition,  high-purity 
nitrogen  was  used  as  carrier  gas.  The  volume  flow  rate  of  nitrogen 
was  80  mL/min.  Non-isothermal  experiment  runs  were  carried  out 
at  heating  rates  of  5,  10,  20,  and  40°C/min  with  temperatures 
ranging  from  ambient  temperature  to  900  °C.  All  samples  were 
heated  to  110  °C  and  held  20  min  (until  a  constant  mass  at  this 
temperature  was  obtained)  to  provide  moisture  removal.  During 
heating,  sample  mass  and  furnace  temperatures  were  recorded. 
Each  experiment  was  repeated  at  least  twice  for  repeatability. 


2.3.  Kinetic  analysis 


The  kinetics  data  obtained  from  TGA  are  crucial  for  the  efficient 
design  and  operation  of  the  pyrolysis  process  (Liang  et  al.,  2014). 
As  a  kind  of  heterogeneous  solid  state  reaction,  general  non- 
isothermal  solid  decomposition  reaction  rate  can  be  expressed  as 
in  Eq.  (1); 

dx 

0) 


where  x  is  conversion  and  given  by; 

W0-W 
X  ~  W0  -  W„ 


(2) 


In  Eq.  (2)  W0,  W  and  refer  to  initial,  instantaneous  and  final 
weights,  respectively.  The  temperature  dependent  reaction  rate 
constant  k(T)  is  expressed  by  Arrhenius  equation  as; 

k(T)=Aexp(-E/RT)  (3) 


where  A  is  the  pre-exponential  factor,  E  is  the  activation  energy  of 
the  reaction,  R  is  the  universal  gas  constant,  and  T  is  the  absolute 
temperature. 

Substituting  Eq.  (3)  into  Eq.  (1)  gives; 

^  =  Ae(-%(x)  (4) 


where  x  is  conversion;  t  is  reaction  time;  A  is  pre-exponential  fac¬ 
tor;  E  is  activation  energy;  R  is  the  universal  gas  constant;  T  is  the 
absolute  temperature.  f(x)  represents  reaction  model  which  is  a 
function  of  x. 

Taking  the  logarithm  of  both  sides  in  Eq.  (4)  yields  the  Friedman 
equation,  Eq.  (5),  where  the  activation  energy  can  be  calculated 
regardless  of  the  decomposition  reaction  model. 

ln  (dt)  =  lnA  +  ln/M  -  ^  (5) 


The  activation  energy  (E)  can  be  obtained  from  the  slope  of  the 
curve  by  plotting  ln (dx/dt)  against  1/T  under  a  given  value  of  x. 

Model  free  Kissenger-Akahira-Sunose  (ICAS)  method  was  also 
applied  on  the  TGA  data  to  obtain  E.  According  to  the  KAS  method, 
E  can  be  determined  from  the  slope  of  the  lines  generated  through 
the  linear  correlation  of  ln  (PIT2)  with  1/T  under  a  given  value  of  x, 
as  shown  in  Eq.  (6)  (Gao  et  al.,  2013): 


ln 


(*)  =  in 

AE 

\T2J 

kwJ 

E 

Rf 


(6) 


where  K  shows  the  heating  rate  (°C/min). 

Model  free  methods  are  iso-conversional  methods  where  the 
activation  energy  is  a  function  of  the  conversion  (Lopez- 
Velazquez  et  al.,  2013).  These  methods  are  more  reliable  than 
model  fitting  methods  in  the  determination  of  the  activation 
energy  because  they  can  be  used  without  knowledge  on  the  reac¬ 
tion  model.  In  this  study,  we  applied  two  kinds  of  model  free  meth¬ 
ods,  Friedman  and  KAS,  to  obtain  the  activation  energy  and  kinetic 
parameters.  The  intention  of  choosing  those  two  types  of  model- 
free  methods  was  to  evaluate  and  verify  the  consistency  of  our 
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results  through  two  different  algorithms:  the  differential  one  of 
Friedman  and  the  integral  one  of  KAS. 

The  Coats-Redfern  method,  which  has  been  widely  applied  to 
investigate  the  thermal  decomposition  was  used  in  this  study  to 
obtain  the  parameters  of  thermal  events.  Mean  E  value  from  model 
free  methods  was  used  in  Coats  Redfern  equations  to  determine 
kinetic  parameters,  the  pre-exponantial  factor  (A)  and  reaction 
order  (n)  (Damartzis  et  al.,  2011;  Kantarelis  et  al.,  2011). 

The  Coats-Redfern  equations  can  be  outlined  as  following  (Gao 
etal.,  2013); 


Table  1 

Proximate  and  ultimate  analysis  and  higher  heating  values  of  PE  and  some  other 
macroalgae. 


PE 

S.  japonicab 

Sargassum  sp.c 

Proximate  analysis  (%) 

Moisture 

11.55 

6.90 

9.34 

Volatile  matters 

48.20 

68.79 

44.50 

Ash 

27.45 

20.21 

36.82 

Fixed  carbon 

12.80 

4.10 

9.34 

Ultimate  analysis 

C 

35.81 

32.89 

26.70 

H 

5.93 

6.17 

4.23 

N 

6.86 

0.93 

1.35 

O 

51.40a 

60.01 

67.53 

Higher  heating  value  (MJ/kg) 

12.54 

12.11 

10.10 

a  By  calculation. 
b  Kim  etal.  (2012). 
c  Kim  et  al.  (2013). 


(a)  TGA 


i-a-*)1-" 

T2(l  -  n 


In 


AR 

M\ 


Ea 

RT 


(for  n 


and, 


In 


ln(l  -  x)' 
T2 


Ea_ 

RT 


(for  n  =  1) 


(12) 


(13) 


where  n,  T,  A,  R,  )S,  and  E  are  the  reaction  order,  absolute  tempera¬ 
ture,  pre-exponential  factor,  gas  constant,  heating  rate,  and  activa¬ 
tion  energy,  respectively.  In  the  plot  of  In  [’  r^-n) "]  versusr  (f°r 
n  #  1 )  or  In  ln(*ra)]  versus  j  (for  n  =  1 )  slope  gives  -£/R.  The  inter¬ 
cept  can  be  arranged  as  In  |^§J  where  A  can  be  calculated  (Damartzis 
etal.,  2011). 


Table  2 

Decomposition  characteristics  of  PE  at  different  heating  rates. 


Heating 
rate  (°C/min) 

T,  (°C) 

Emax  (°C) 

Tf(°  C) 

DTG  max 

Mass  loss'3  (%) 

5 

226.20 

262.26 

371.70 

0.21 

49.72 

10 

228.12 

271.54 

402.73 

0.41 

50.22 

20 

228.01 

284.33 

436.46 

0.78 

50.45 

40 

227.91 

301.54 

484.15 

1.81 

51.03 

DTCmax  -  maximum  mass  loss  rate. 

Tj  -  initial  decomposition  temperature  of  main  decomposition. 
Tmax  -  temperature  at  DTG  max- 

Tf  -  initial  decomposition  temperature  of  main  decomposition. 
a  Mass  loss  at  main  degradation  stage. 


Fig.  1.  (a)  TGA  curves  of  PE  at  different  heating  rates  and  (b)  DTG  curves  of  PE  at  Fig.  2.  Plots  for  determination  of  activated  energy  at  different  conversion  by 
different  heating  rates.  (a)  Friedman  method  and  (b)  KAS  method. 
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3.  Results  and  discussion 

3.J.  Characterisation  ofP.  elongata  (PE) 

The  results  of  the  proximate  analysis,  ultimate  analysis  and 
higher  heating  value  of  PE  and  the  comparison  with  some  macro- 
algaes  previously  studied  can  be  seen  in  Table  1.  The  characteris¬ 
tics  of  PE  were  similar  with  other  macroalgae.  The  volatile 
matter  content  of  PE  was  lower  than  that  of  S.japonica  but  slightly 
higher  than  that  of  Sargassum  sp.  The  ash  content  was  27.45%  for 
PE  and  it  was  between  20.21%  and  36.82%  for  S.  japonica  and 
Sargassum  sp.,  respectively.  The  low  ash  content  is  important 
because  high  ash  content  can  cause  aggregation  in  processes  and 
yield  inefficient  heat  transfer  rates.  It  can  be  seen  from  results  that 
PE  is  a  carbon  and  oxygen  rich  biomass  containing  5.93%  H  and 
6.86%  N.  The  empirical  formula  of  PE  can  be  expressed  as 
CH1.99Ch.08No.i64-  The  high  heating  value  (HHV)  of  PE  (on  dry 


biomass)  was  measured  as  12.54  MJ/kg  which  is  slightly  higher 
than  that  of  S.  japonica  and  Sargassum  sp.  The  higher  carbon  and 
hydrogen  content  can  be  the  reason  for  the  higher  value  of  HHV. 

3.2.  Thermogravimetric  analysis 

The  mass  loss  (TG)  and  differential  mass  loss  (DTG)  thermo¬ 
grams  of  the  thermal  decomposition  of  PE,  at  four  heating  rates, 
5,  10,  20  and  40  °C  min-1  under  nitrogen  atmosphere,  are  shown 
in  Fig.  la  and  b.  When  the  thermogravimetric  analysis  was  con¬ 
ducted,  the  final  weights  were  49-51%  for  each  heating  rate  at 
650  °C.  Thus,  the  char  content  is  expected  to  be  approximately 
50%  indicating  the  existence  of  a  large  amount  of  inorganic  com¬ 
pounds  (Kim  et  al.,  2013).  The  main  weight  loss  occurred  in  the 
range  of  225-400  °C,  with  a  main  peak  in  the  DTG  at  280-300  °C 
and  a  shoulder  approximately  between  300-400  °C.  Such  trends 
may  be  attributed  to  the  decomposition  of  carbohydrates 


Table  3 

The  fitted  equations,  correlation  coefficients  (R2)  and  activation  energies  (£,  kj/mol)  obtained  by  the  Friedman  and  KAS  methods. 


Conversion  (x) 

Friedman  method 

KAS  method 

Fitted  equation 

R2 

£ 

Fitted  equation 

R2 

Ea 

0.10 

y  =  -17550x  + 27.757 

0.9974 

147.57 

y- 

-15432X  + 18.987 

0.9823 

128.30 

0.15 

y  =  -18253x  + 28.78 

0.9967 

151.76 

y= 

-15799X  + 19.132 

0.9856 

131.35 

0.20 

y  =  -19025x  + 29.835 

0.9901 

158.17 

y- 

-16168X  + 19.381 

0.9803 

136.91 

0.25 

y=-19984x  + 31.131 

0.9901 

166.15 

y= 

-16995X  + 20.492 

0.9891 

141.30 

0.30 

y=-21090x  + 32.475 

0.9805 

175.34 

y- 

-17490X  + 20.946 

0.9841 

145.41 

0.35 

y  =  -21204x  + 31.92 

0.9747 

176.29 

y- 

-18069X  + 21.448 

0.9695 

150.23 

0.40 

y  =  -18581x  + 26.564 

0.9899 

154.48 

y= 

-17203X+ 19.325 

0.9707 

143.03 

0.45 

y  =  -16608x  + 22.483 

0.9997 

138.08 

y= 

-16935X  + 18.255 

0.9809 

140.80 

0.50 

y  =  _15144x  + 19.291 

0.9966 

125.91 

y- 

-15893X+  15.728 

0.9979 

132.13 

0.55 

y=-12846x  + 14.807 

0.9522 

103.81 

y= 

-14364X  + 12.685 

0.9875 

119.42 

0.60 

y  =  -11121x  + 11.16 

0.9514 

92.46 

y- 

-13071X  + 9.8923 

0.9728 

108.67 

0.65 

y  =  -10063x  + 6.0171 

0.9980 

83.66 

y= 

-9975.5X  + 4.1693 

0.9900 

82.94 

0.70 

y  =  -9662.1x  + 7.7476 

0.9993 

80.33 

y- 

-9030.9X  + 1.9184 

0.9949 

75.08 

0.75 

y  =  -8667.7X  +  4.4408 

0.9894 

72.06 

y- 

-6837.6X-  1.7174 

0.9928 

56.85 

0.80 

y  =  -8533.3x  + 5.0813 

0.9757 

71.11 

y- 

-6144.1X- 3.1219 

0.9903 

51.08 

5°C/min 


10°C/min 


40°C/min 


Fig.  3.  Coats-Redfern  plots  for  PE  at  different  heating  rates. 
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(225-300  °C)  and  proteins  (300-400  °C)  (Bae  et  al.,  2011).  Similar 
results  were  also  observed  during  the  pyrolysis  of  some  other  mac¬ 
roalgae  (Bae  et  al.,  2011;  Kim  et  al„  2012,  2013;  Li  et  al.,  2011). 
Additionally,  as  seen  in  Fig.  la,  above  500  °C  a  slight  weight  loss 
occurred  due  to  further  devolatilization  of  the  formed  biochar.  This 
phenomenon  can  be  the  result  of  the  breakdown  of  C-C  and  C-H 
bonds  of  the  biochar  (Liang  et  al.,  2014).  As  seen  in  Fig.  la  and  b, 
both  TGA  and  DTG  curves  showed  similar  trends  for  all  heating 
rates.  However  a  slight  shift  with  the  increase  in  heating  rate 
was  observed.  The  heating  rate  effected  the  pyrolysis  rate  and 
decomposition  temperatures  during  pyrolysis.  It  also  appears  that 
as  the  heating  rate  was  increased,  the  thermal  degradation  process 
was  delayed.  At  higher  heating  rates,  the  samples  reach  the  tem¬ 
perature  at  a  shorter  time  as  the  decomposition  temperatures  were 
shifted  to  higher  values  (Wu  et  al.,  2014).  The  heating  of  biomass 
particles  occurs  more  gradually  at  lower  heating  rates,  thus 
enabling  an  improved  and  more  effective  heat  transfer  to  the  inner 
portions  and  among  the  particles  (Idris  et  al.,  2012;  Kim  et  al., 
2013;  Gai  et  al.,  2013;  Zhao  et  al.,  2013).  Therefore,  at  higher 
heating  rates,  the  decomposition  process  proceeded  slower  due 
to  inefficient  heat  transfer. 

Table  2  represents  the  decomposition  characteristics  of  PE  pyro¬ 
lysis  reaction  at  different  heating  rates.  As  seen  in  Table  2,  the  main 
pyrolysis  process  proceeds  in  a  range  from  about  225  °C  to  372  °C 
for  low  heating  rates  and  484  °C  for  high  heating  rates  leading  to 
formation  of  volatile  matter  and  the  biochar.  It  can  be  seen  both 
from  Table  2  and  Fig.  la  and  b  that  the  maximum  mass  loss  rate 


Table  4 

Kinetic  parameters  for  PE  biomass  pyrolysis  in  nitrogen  atmosphere. 


Heating  rate  (°C/min) 

n 

A 

R2 

5 

7.6 

4.88  x  10” 

0.9914 

10 

8.1 

1.07  x  1012 

0.9908 

20 

8.4 

4.47  x  10” 

0.9872 

40 

7.4 

5.00  x  10" 

0.9854 

(a)  5°C/min 


T  (°C) 


(c)  20°C/min 


T  (°C) 


(DTGmax)  increases  with  the  increment  in  the  heating  rate.  When 
heating  rate  was  5  °C/min,  the  DTGmax  was  0.21  mg/min.  However, 
DTGmax  increased  to  1 .81  mg/min  when  the  heating  rate  was  40  °C/ 
min.  The  DTGmax  is  a  function  of  the  heating  rate  and  increases 
with  the  increase  in  heating  rate  due  to  a  higher  amount  of  thermal 
energy  which  promoted  the  heat  transfer  between  the  surround¬ 
ings  and  interiors  of  the  sample  (Kim  et  al.,  2013).  The  peak  tem¬ 
peratures  were  also  influenced  by  the  heating  rate.  As  the 
heating  rate  increased  the  maximum  peak  temperature  was 
slightly  shifted  to  higher  temperatures.  This  can  be  the  cause  of 
heterogeneous  structure  of  macroalgae.  At  high  heating  rates, 
decomposition  of  carbohydrate  and  protein  components  of  macro¬ 
algae  may  occur  simultaneously.  As  seen  in  Fig.  2b,  at  a  heating 
rate  of  40  °C/min,  shoulder  peak  disappeared  and  an  overlapping 
peak  was  observed  for  decomposition  of  PE  (Liang  et  al.,  2014). 

3.3.  Kinetic  analysis 

The  model  free  Friedman  and  KAS  methods  were  used  to  eval¬ 
uate  the  apparent  activation  energy.  Methods  were  employed  at 
different  heating  rates  ranging  from  5  to  40  °C/min,  for  fractional 
conversions  varying  from  0.10  to  0.80  to  determine  the  variation 
of  the  apparent  activation  energy  during  the  thermal  decomposi¬ 
tion  process.  The  regression  lines  obtained  using  Friedman  and 
KAS  methods  are  shown  in  Fig.  2a  and  b,  respectively.  The  fitted 
equations,  correlation  coefficients  ( R 2)  and  activation  energies  are 
summarized  in  Table  3.  R2  values  were  higher  than  0.9700  for  all 
lines.  Therefore,  it  can  be  concluded  from  Fig.  2  and  Table  3  that 
the  model  free  methods  used  are  reliable  in  determining  the  acti¬ 
vation  energy.  The  activation  energy,  E,  based  on  Eq.  (5)  can  be 
obtained  from  a  plot  of  ln(dx/dt)  vs.  1/T.  Fig.  2a  presents  plots  of 
ln(dx/dt)  vs.  1/T  at  various  conversions  ranging  from  0.10  to  0.80. 
Similarly,  as  seen  in  Fig.  2b,  plotting  the  left  side  of  Eq.  (6), 
In  vs.  inverse  temperature,  1  /T,  for  different  heating  rates, 

straight  lines  can  be  produced.  As  seen  in  Fig.  2b,  the  activation 
energy  can  be  estimated  by  the  slope  of  produced  straight  lines. 


(b)  10°C/min 


T  (°C) 


(d)  40°C/min 


T(°C) 


Fig.  4.  Simulation  of  PE  pyrolysis  using  the  kinetic  data  calculated  for:  (a)  5,  (b)  10,  (c)  20  and  (d)  40  °C/min. 
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Table  5 

Comparision  of  kinetic  parameters  obtained  from  pyrolysis  reactions  of  various 
marine  biomass  residues. 


Marine  biomass 

E  (kj/mol) 

A  (s  ’) 

n 

Reference 

C.  vulgaris 

332.6 

1.3  x  1028 

9.0 

Chen  et  al.  (2012) 

S.  platensis 

97.0 

9.16  x  107 

1.98 

Peng  et  al.  (2001) 

C.  protothecoides 

52.5 

1.11  x  104 

1.88 

Peng  et  al.  (2001) 

D.  tertiolecta 

146.067 

2.28  x  1013 

2.4 

Shuping  et  al.  (2010) 

D.  tertiolecta 

171.85 

1.09  x  1014 

5.3 

Wu  et  al.  (2014) 

E.  prolifera 

178.89 

6.321  x  10'° 

5.4 

Wu  et  al.  (2014) 

P.  elongata 

121.355a 

6.263  x  10"3 

7.9a 

This  study 

3  Mean  values. 


The  activation  energies  calculated  by  KAS  and  Friedman  meth¬ 
ods  were  slightly  different.  However,  the  variations  in  the  activa¬ 
tion  energy  with  proceeding  conversions  were  similar  in  both 
methods.  The  activation  energy  was  increased  up  to  0.35  conver¬ 
sion  then  it  was  decreased.  The  activation  energy  was  increased 
from  147.57  to  176.29  kj/mol  and  decreased  to  71.11  kj/mol  in 
the  Friedman  method.  It  increased  from  128.30  to  150.23  kj/mol 
and  then  decreased  to  51.08  kj/mol  in  the  KAS  method.  This  trend 
can  be  caused  by  the  degradation  of  carbohydrate  and  protein  con¬ 
tent  of  macroalgae  PA  (Slopiecka  et  al.,  2012).  Mean  values  were 
126.48  kj/mol  and  116.23  kj/mol  for  Friedman  and  KAS  methods, 
respectively.  The  mean  of  these  values  were  accepted  as  the  mean 
activation  energy  and  used  in  Eqs.  (12)  and  (13)  for  the  calculation 
of  kinetic  parameters,  including  pre-exponential  factor  and  reac¬ 
tion  order.  As  seen  from  Fig.  3,  for  all  heating  rates,  R2  values  for 
Coats-Redfern  plots  were  higher  than  0.9800. 

The  obtained  pre-exponential  factor  and  reaction  order  values 
at  different  heating  rates  are  shown  in  Table  4.  It  should  be  noted 
that  these  reaction  orders  which  were  obtained  by  iso-conversion- 
al  methods  have  no  physical  meaning  and  can  be  considered  as  fit¬ 
ting  parameter  (Damartzis  et  al.,  2011;  Kantarelis  et  al.,  2011; 
Liang  et  al.,  2014).  The  obtained  kinetic  parameters  were  used  to 
simulate  the  pyrolysis  process  of  alga  PE.  As  seen  from  Fig.  4, 
model  and  experimental  data  were  in  a  good  agreement. 

A  comparison  of  the  kinetic  parameters  of  pyrolysis  for  PE  and 
some  other  marine  biomasses  are  shown  in  Table  5.  The  activation 
energy  for  pyrolysis  of  PE  is  similar  to  the  reported  values  for 
D.  tertiolecta,  E.  prolifra.  However,  the  activation  energy  for 
C.  vulgaris  was  higher.  In  contrast,  the  activation  energies  for  S. 
plantensis  and  C.  protothecoides  were  lower  than  that  of  PE.  An 
obvious  difference  also  can  be  seen  in  the  pre-exponential  factors 
and  reaction  orders.  The  differences  in  the  kinetic  parameters 
may  be  due  to  the  differences  in  the  composition  of  the  biomasses. 
Also,  the  usage  of  different  experimental  conditions  and  various 
methods  for  calculations  may  cause  that  the  derived  kinetic 
parameters  to  differ  from  each  other  (Slopiecka  et  al.,  2012; 
Liang  et  al.,  2014). 

4.  Conclusion 

The  incorporation  of  marine  biomass,  PE,  into  energy  produc¬ 
tion  depends  on  understanding  of  the  kinetics  of  pyrolysis  and 
degradation  reactions.  As  a  readily  available  and  abundant  raw- 
material,  activation  energies  for  PE  were  comparable  to  those  of 
previously  studied  algae  residues.  A  variance  in  the  activation 
energy  with  proceeding  conversions  was  observed  which  revealed 
that  the  pyrolysis  of  PE  was  composed  of  multi-step  kinetics.  This 
finding  was  consistent  with  the  reported  studies  on  some  other 


algae  species,  D.  tertiolecta,  E.  prolifra,  S.  plantensis  and  C.  prototh¬ 
ecoides.  A  reaction  model  was  established  using  obtained  kinetic 
parameters  and  results  were  in  a  good  agreement  with  experimen¬ 
tal  data. 
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